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A systematic investigation is presented of grain boundaries and 
grain boundary networks in two dimensional flexible membranes with 
crystalline order. An isolated grain boundary undergoes a buckling 
transition at a critical value of the elastic constants, but, contrary to 
previous findings, the shape of the buckled membrane is shown to be 
asymptotically flat. This is in general not true in the case of a network 
of intersecting grain boundaries, where each intersection behaves as the 
source of a long range stress fleld. Unlike the case of isolated disloca- 
tions or disclinations, the energy associated with these stresses is flnite; 
they can, however, destabilize the flat phase. The buckled phase can 
be modeled by an Ising spin-glass with long range antiferromagnetic in- 
teractions. These flndings may be relevant to the understanding of the 
wrinkling transition in partially polymerized vesicles reported in recent 
experiments. 

PACS numbers: 87.22. Bt, 46.30.Lx 
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I. INTRODUCTION 

Flexible membranes are two dimensional (2D) generalizations of linear polymer chains. 
The properties of a 2D membrane, embedded in three dimensional space, depend strongly 
on the internal order, crystalline, hexatic, or fluid. As in other realizations of 2D matter, 
defects, and their interactions, affect crucially the stability of a given phase. For example, 
the melting transition in two dimensions can be described as proliferation of topological 
defects [1] . A membrane, however, is not confined to a plane. Thus, although the stable 
phase of a defect-free crystalline membrane at low temperature is flat, strains induced 
by a defect, such as a dislocation, can be accommodated by displacements in the normal 
direction, resulting in the buckling of the membrane [2] . This process entails a trade-off of 
in-plane stretching energy for curvature energy, and, hence, it occurs when 

-—>^. (1) 

Here, Kq is Young's modulus, n is the bending rigidity, I is a length scale, and 7 is a 
dimensionless constant of order 10^. In membranes of size R, I — R for disclinations, and 
I = y/Rb for a dislocation with Burger's vector b. Thus, these defects always buckle in 
sufficiently large membranes, irrespective of the value of the elastic constants. 

This conclusion, however, does not hold for finite energy defects such as vacancies, 
interstitials, or tightly bound dislocation pairs. In this case, I is of order a lattice constant, 
and thus, the stability of the fiat phase in the presence of such defects is determined by 
the actual value of the elastic constants. This brings about the interesting possibility of a 
buckling transition in an infinite system as a function of temperature. Because k typically 
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increases with temperature, while Kq decreases, the condition (1) may be associated with 
a boundary in the (/^, r)-plane, separating buckled and flat regimes (see Fig. 1). 

Recent studies of membranes with defects and quenched random disorder have been 
prompted in part by experiments of Mutz et al. [3], showing that partially polymerized 
vesicles undergo a (possibly first-order) reversible phase transition from a high temperature 
phase characterized by a smooth, "floppy" surface, to a low temperature phase where the 
vesicles appear rigid and highly wrinkled. This wrinkling transition has been likened to 
the spin-glass transition of magnetic systems [3,4]. To understand the mechanism of the 
transition, it is helpful to observe that the effect of the partial polymerization is to nucle- 
ate microcrystalline domains around polymerized patches. Upon cooling, crystalline order 
arises in the unpolymerized material surrounding these domains. This has led to specu- 
lation that the grain boundaries, resulting when microcrystals with different orientations 
meet, may be responsible for the observed wrinkled structure. The following argument 
was given to support this scenario [5] . 

Consider a low angle grain boundary in a 2D crystal. This can be viewed as a row of 

— * 

dislocations, with Burger's vector h perpendicular to the boundary. Let 6 be the boundary 
tilt angle. The spacing between dislocations, /i, is then given by Frank's law [6], h = 
6/2 sin(^/2) h/9. On length scales large compared to /i, the stress flelds of the dislocations 
are screened exponentially, provided the average Burger's vector is strictly perpendicular; 
on shorter length scales, however, the behavior of the grain boundary should be dominated 
by the individual dislocations. Computer simulations of tethered membranes [2] have 
shown that isolated dislocations in a finite membrane of size R destabilize the fiat phase, 
provided that R > 127K,/Kob = Rb, consistent with Eq. (1) with I = \fRb. Hence, one 

3 



might expect that, if h > Rb, the grain boundary will buckle out of flat space, possibly- 
taking the shape of a rippled hinge, possibly with large out-of-plane displacements, which 
would then appear as wrinkles. 

The present work is devoted to a systematic study of grain boundaries and grain 
boundary networks in two dimensional flexible membranes with crystalline order. While 
we flnd that an inflnite grain boundary does indeed undergo an intriguing second order 
transition from a flat to a corrugated phase as the ratio of the bending rigidity to the 
shear modulus or tilt angle is varied, we show that the out-of-plane displacement decays 
exponentially away from the boundary, i.e., the shape of the membrane remains asymp- 
totically flat. We also flnd, however, that this conclusion does not hold in the case of a 
network of intersecting grain boundaries, where each intersection behaves as the source 
of a long range stress held, e.g., as a disclination multipole. Unlike the case of isolated 
dislocations or disclinations, the energy associated with these stresses is generally flnite; 
they can, however, cause buckling at a sufficiently low temperature, determined by the size 
I of the (localized) multipole source. 

Thus, upon cooling of a polycrystalline membrane, in which a grain boundary network 
has formed, one expects a hierarchy of buckling transitions to take place, with the low- 
angle grain boundaries puckering first, followed by the high-angle ones (low-angle grain 
boundaries correspond to very large values of P in Fig. 1). Eventually, as the network of 
buckled grain boundaries begins to percolate across the membrane, the (point-like) grain 
boundary nodes will buckle, causing large wrinkles to develop fully. We will show that the 
long range, often highly anisotropic, interactions in the buckled phase can favor buckling 
of a pair of nodes in opposite directions. Buckling of these nodes can be described by 
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randomly placed Ising-like degrees of freedom, with antiferromagnetic interaction. Thus, 
frustration will play an important role in the conformation of these membranes in the 
buckled phase, much as in spin systems with (long range) random antiferromagnetic ex- 
change [7,8]. Ising spin-glass behavior seems likely at low temperature, which may provide 
a plausible explanation of the wrinkling transition reported in Ref. [3] . 

The balance of our paper is as follows. We review some results of the theory of 
membrane buckling in Section II, and establish a connection with the Landau theory of 
second order transitions. In Section III, we study isolated grain boundaries and grain 
boundary networks within the continuum elastic theory of membranes. Our findings are 
supported by numerical simulations presented in Section IV, and we conclude in Section 
V with a discussion of simplified models based on this work. 



We begin by reviewing the continuum theory of elastic membranes [9] . Let the shape 
of a membrane be parameterized by a three dimensional displacement vector with respect 
to a reference plane {x,y), and let u, f be the in- plane and out-of- plane components 
respectively. The energy of the membrane in a given configuration is the sum of a stretching 
and a bending term. 



Summations over repeated indices are understood. Here k, is the bending rigidity, A and /i 



II. BUCKLING OF ELASTIC MEMBRANES 



F = F, + Ft, 




(2) 
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Lame coefficients, and the strain tensor is given by 

1 



u 



= 7^(^i%- + djUi + difdjf). (3) 



Note the manifest up-down symmetry (/ — > — /). A given configuration of a fluctuating 
membrane occurs with probabihty proportional to exp(— F//cbT). 

Minimizing the energy with respect to the displacement yields the equations 

nV^f = aijdidjf 

(4) 

diaij = 0, 

where the stress tensor cry is related to the (nonlinear) strain tensor by 

aij = 2iiUij + XuiiSij. (5) 

Because is divergenceless and symmetric, it can be expressed through a single scalar 
potential, the Airy stress function x- 

d^X _ d'^X _ d^X 

Equations of motion for the energy functional are most conveniently derived in terms 
of / and X- Cartesian coordinates, the equations, obtained originally by von Karman, 
read 

«VV ^ + _ (7a) 

dy'^ dx"^ dx^ dy"^ dxdy dxdy 

Here Kq = 4)Li(/i + A)/ {2^ + A) is the two-dimensional Young's modulus, and the source 
term, 5'(r), is the total density of defects. For an isolated disclination at tq with charge s, 
one has 

S{f^ = s5{r-ro), (8a) 
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while for a disfocation with Burger's vector b we have 

S{r)^e,jhdj6{f-fo). (86) 
An isolated vacancy, interstitial, or impurity atom [4] is associated with a source term 

where a = X/ (2// + A) is the 2D Poisson ratio and Qq the area change associated with the 
defect. Other defects may be regarded as dislocation or disclination multipoles. 

An effective free energy functional for the normal displacement / can be obtained by 
integrating out the in-plane displacement field u. The result is [10] 



(9) 



where K{r) = {d'^f/dx^){d'^f/dy'^) — {d"^ f / dxdyf is the gaussian curvature of the mem- 
brane at point r. 

This clearly displays that the free energy is the sum of two competing terms. One, 
proportional to Young's modulus Kq, favors configurations where the membrane buckles 
out of the plane to screen out the source term. The other, proportional to the bending 
rigidity, favors configurations with minimum mean curvature, that is, the flat state. A 
buckling instability is then expected to occur, for a given source term, as the ratio Kq/ n 
increases, or, for given elastic constants, as the strength of the source increases. Because 
the stretching energy of isolated dislocations and disclinations in flat space diverges with 
system size, buckling occurs immediately for these defects for all values of the elastic 
constants in an infinite system. 



The von Karman equations are extremely diflBcult to solve, and most of what we 
know on the buckling of isolated defects was obtained by numerical simulations. One can, 
however, make rigorous arguments about the buckling transition, mainly through linear 
stability analysis. Observe that, if the energy of a given configuration S{r) is finite in flat 
space, then for small but finite Kq the fiat phase will be stable. Then / = 0, and one can 
solve for the Airy stress function of the fiat state Xfi 



Stability with respect to out-of-plane buckling then demands that the linear eigenvalue 
problem 



admit only nonnegative eigenvalues rj. Equivalently, we see that, as the ratio Kq/k in- 
creases, a buckling instability will set in as soon as an eigenvalue rj becomes negative 
(?7 = 0~). To this eigenvalue corresponds a stable configuration /^(r) ^ 0: the up-down 
symmetry is now broken. Thus, fn{r) can be regarded as the order parameter of a Landau 
theory that describes the buckling transition. The Landau free energy, given by Eq. (9), 
contains terms quadratic and quartic in /. Suppose for simplicity that Kq is fixed (it 
can always be scaled out of the problem). Then the buckling instability is reached by 
decreasing k from infinity to a critical value, k*. 

It is easy to see that the behavior of f{r) near the critical point is that expected 
from Landau theory. Let fn{r) be the normalized eigenfunction of Eq. (11) obtained 
at the critical point n = n* . Below the buckling transition, but sufficiently close to it 



V\f=KoS{r). 



(10) 




(11) 
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(k = K* — 5k,, 5k, <^ k,*), we look for a solution of the von Karman equations of the form 

Xir) = Xfir) + ^Xir) 

(12) 

/(r) = afnir), 

where a is a scalar Ising-like order parameter. We immediately find 

5x{r) = -a^KoV-^k{r), (13) 

where K{r) is the Gaussian curvature of fr){r)- Multiplying Eq. (11) on the left hy fn{r) 
and integrating over all space, we find (putting r(r) = i^oV~^K(r)) 

and, thus, a oc (5k) 

III. GRAIN BOUNDARY BUCKLING 

Grain boundaries are interfaces joining continuously two crystalline domains, A and B, 
of different orientation (see Fig. 2). Consider a reference lattice, e.g., the infinite extension 
of domain A. Domain B can be obtained by rotating A through an angle 9 about the 
vertical. Consider a closed contour in the reference lattice. If the part that lies in B is 
rotated as described above, and the contour encloses a segment of length h of the grain 
boundary, the contour will fail to close by an amount [6] 

b = 2hsm-. (15) 

Thus, dislocations must lie along the grain boundary, with Burger's vector density deter- 
mined by the tilt angle 9. Dislocations need not all have the same Burger's vector (so 
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long as it has the proper magnitude and direction on average), but in the simplest case 
of a symmetric grain boundary, the boundary is equivalent to an infinite row of identical, 
equally spaced dislocations. We begin by studying this case in detail. 

Consider a membrane with a density of defects consisting of a row of dislocations, 
lying on the |/-axis, with Burger's vectors b: 

(ft ft \ °^ 

E ^(y-rih). (16) 

The Airy stress function in flat space is obtained from Eq. (10), 

II 1 °° 1 

X(0 = ~^y^^^f^ + ^^^°4^ 5-/ n ^^P(~27rn|a;|//i)a;cos(27rn?///i) 

n=l 

1 °° 1 f h \ 

—bxKo— y exp(— 27rn|a;|//i) I - — h n\x\ I sm{27rny / h) . 
An ^—^ n V 27r / 

n=l ^ ^ 

The requirement that the energy per unit length of the boundary be finite implies, in 
fiat space, the well known condition that there be no net component of the Burger's vector 
parallel to the boundary; i.e., by = 0. The stresses then decay exponentially in the grains 



(17) 



as \x\ — > oo, 



Kob f 271 \f 27r A . f2n 
exp — — 1 + — sm —y 



2h "V ^ 7 V 7 V ^ 

Kob f 27r, A / 27r, A . f 2n , 
"^yy - ^exp —r\A 1 - -r\x\ sm —y (18) 



2h \ h ' 'J \ h ' ' J \ h 

Kob ( 27r, A (2ti 

CT-Eu ~ — — exp \x\ he cos — i 

^ 2h ^ 7 V ^ ' 

This expression for the stress tensor can be used to perform a stability analysis according 



to Eq. (11). The buckling problem becomes simpler, and more physically intuitive, if we 
work in a zero-range approximation, which should be particularly appropriate for high- 
angle grain boundaries, where h is of order a few lattice constants. We expect that this 
approximation describes the asymptotic physics more generally in the limit x ^ h. In this 
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approximation, which preserves the integrated stress, the components of the stress tensor 
become 



^yy ^ '-^xy ~ 0- 



(19) 



Figure 3 illustrates the physical origin of the stresses above, for the case of a 21.8° sym- 
metric grain boundary in a triangular lattice. 

Using Eqs. (6) and (19), we see that the stability problem now requires finding non- 
trivial solutions of the equation 

V/-^5(x)sm(-y)^ = 0. (20) 

We look for a periodic solution in the direction of the boundary, with period h. Note that 
/ must be of the form (up to an additive constant and a linear term in x) 

fi'^ = ~f^ ^ /(?n)exp(zg^y)exp(-|Q„a;|)(l + |g^a;|), (21) 

n=—oo 

with Qn = 27rn/h, which shows that, upon buckling, the surface remains asymptotically 
flat as \x\ — > oo. Substituting this expression for / into Eq. (20) and taking the Fourier 
transform, one obtains the following recursion relation for q^fiQu) = fiQn)' 

f{qn+i) + -Tr^\qn\f{qn) + f{qn-i) = 0. (22) 

A buckling instability occurs when the determinant of the tridiagonal matrix defined by 
the recurrence first vanishes as the amplitude of the diagonal terms (i.e., of k) is decreased 
from infinity. This happens for 
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which yields the estimate of the critical value of the elastic constants for the buckling of 
the zero-range model of symmetric grain boundaries 

Kobh 



K 



190. (24) 



For fixed elastic constants, this can be regarded as a condition on the tilt angle, 9 
b/h. Just below the flat space Kosterlitz-Thouless melting temperature we have Kob'^ = 
IQttTkt, so grain boundaries will buckle for T < Tkt provided 9 < 9c ^ Q.27Tkt/ 

The same conclusion, concerning the shape of the membrane after buckling, could be 
arrived at qualitatively from Eq. (2). Indeed, we note that, if buckling is to reduce the 
stretching energy, then the gradient of / ought to approximately cancel the components 
of the in- plane displacement [see Eq. (3)]. Then and, thus, / itself, should decay 

exponentially (i.e., as the square root of the in-plane stress). 

This finding appears in conflict with evidence of a long range displacement / reported 
in Ref. [5] , where a paper model of a grain boundary in a triangular lattice was constructed. 
The model is pertinent to the limit Kobh/n » 1 because of the large in-plane shear 
modulus. The grain boundary buckles, forms a hinge-like structure with a definite dihedral 
angle, and displays intriguing ripples in one of the two grains. Here, we show that finite size 
effects and the slight misalignment of the Burger's vectors with the grain boundary normal 
affect the conclusions crucially. Figure 4 displays an atomistic view of the grain boundary 
studied in Ref. [5] . Note that one domain is tilted by the same angle as each of the grains 
of Fig. 3, but the other is not. Therefore, it is no longer possible to create a stress-free 
configuration with coincident atomic sites: to bring the grains to match at the boundary 
in flat space, one has to compress uniformly (i.e., for all x) the untilted grain by a factor 
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-\/7/3 in the y direction (or to stretch the tilted grain correspondingly). If the grain is 
allowed out of the plane, then it can develop ripples in such a way that now the projection 
of the atomic rows on the plane is compressed, but the equilibrium interatomic distance 
in 3D space is restored. In dislocation language, this grain is characterized by Burger's 
vectors with a component parallel to the grain. According to Eq. (17), we have now 
ayy ~ —hyKox / 4h\x\, which is indeed a constant compressive load on the untilted grain, 
and tensile load on the tilted one. When the in-plane stresses are relieved by formation 
of ripples in the compressed grain, stretching energy is traded for bending energy. This is 
certainly advantageous in a paper model, which can be easily bent and hardly stretched; 
if the grains are infinite, however, and if ^ 0, the ripples will themselves cost infinite 
bending energy per unit boundary length. This "infinite" energy can nevertheless still be 
less than the "infinite" energy of the configuration in fiat space, which drives the buckling 
in the finite model. 

Until now, we have dealt with infinite, periodic arrays of dislocations. We emphasize 
that periodicity played a crucial role in our arguments that grain boundary buckling does 
not generate any long range out-of-plane displacement. This may not be, however, the 
most physically relevant situation. In fact, while grain boundaries cannot terminate inside 
a crystalline membrane, three or more of them can radiate outward from a common vertex. 
The breakdown of the periodicity of the arrays has extremely interesting consequences. We 
illustrate them by explicitly considering the case of three semi-infinite arrays of dislocations 
radiating outwards from the origin. 

Let us begin by computing the shear stress produced by a half line of dislocations 
lying on the negative y axis, equally spaced, and with Burger's vector pointing in the x 

13 



direction: 



n=-oo + (y - n/l)2) 

With the aid of the Euler-Maclaurin summation formula [11], 

Y^F{n) = -\{F{N)^F{^))^ dxF{x)+ d J - ^ ^^^^ f (26) 

n=l -^0 -^0 ^ n=l ^ 

we find a multipole expansion for the shear stress field (neglecting exponentially small 
terms) : 



— * 

When grain boundaries meet at a vertex, the Burger's vectors hi and dislocation spac- 
ings hi cannot be arbitrary. This is because they determine the angles, through which the 
crystal orientation is rotated at each boundary (see Frank's law). Imposing the constraint 
that the sum of the Burger's vector densities radiating from the vertex vanish is equiv- 
alent to requiring that no net disclinicity be present at the vertex, i.e., that the sum of 
the monopole terms of the stress field vanish. We have an additional degree of freedom 
at our disposal in the construction of the vertex, which is a rigid shift of the dislocations 
along one arm. This does not affect the monopole moment, but can be used to make the 
unphysical dipole terms (corresponding to a Burger's vector which is not a lattice vector) 
vanish. Higher order multipoles, however, are still present. Below some critical value of 
the elastic constants, the vertex is stable in the flat phase, because the long-range stresses 
it generates are square integrable, and thus, the stretching energy is finite. Above that 
critical value, a buckling transition will take place to a phase with long range out-of-plane 
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displacements. A grain boundary network can thus be regarded as a system of point-like 
objects, the vertices where the grains intersect, with long range, anisotropic interactions. 
The nature of the buckled phase of such a system is highly nontrivial, and numerical 
simulations suggest that spin-glass-like frustration will be important at low temperature. 

IV. NUMERICAL SIMULATIONS 

We have carried out numerical simulations of grain boundary buckling, following the 
method described in Ref. [2] . Briefly, we work with a 2D triangular lattice of atoms, which 
is allowed to bend. In the ideal, unstrained crystal, each atom is sixfold coordinated, and 
each bond has unit length. To represent quenched defects, note that a positive (nega- 
tive) disclination corresponds to a fivefold (sevenfold) coordinated atom, a dislocation is 
a disclination dipole, and so on. The stretching energy is defined as a sum over nearest 
neighbor atoms, 

F, = ^Ko J2{\R,-Rj\-lf, (28) 

<i,j> 

— * 

where R is the 3D atomic coordinate. Each elementary triangle is assigned a unit normal, 
fioi, and the bending energy is defined as a sum over nearest neighbor triangles, 

Fs = ^K J2 {i-na-nf3)- (29) 

<a,f3> 

The total energy is minimized, using a conjugate gradient search [12]. 

First, we have studied the symmetric, high angle grain boundary, whose unrelaxed 
configuration in flat space is shown in Fig. 3. Periodic boundary conditions are im- 
posed along the grain, and free boundary conditions in the transverse direction. Figure 5a 
shows a relaxed configuration after buckling, obtained with a ratio of the elastic constants 
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Kobh/n 300. We have determined numerically that the critical value of the elastic 
constants for this symmetric grain boundary is 

^ ^ 120. (30) 

This value is smaller than predicted by the zero range model, and is very close to that which 
would be obtained for an isolated dislocation in a membrane of finite size h. Note that the 
vertical displacement never exceeds one lattice unit, and that the asymptotically flat state 
is reached very quickly as we move away from the grain boundary. For comparison. Fig. 5b 
shows the results of the relaxation of a finite segment of the same grain boundary; free 
boundary conditions are imposed in all directions. Now one can clearly see a long range 
structure in the out-of-plane displacement, due to finite size effects; the long range stresses, 
caused by terminating the boundary freely, can be estimated from Eq. (27). Figure 5c 
shows the relaxation and buckling of the asymmetric boundary shown in Fig. 4, again with 
periodic boundary conditions along the boundary and free boundary conditions elsewhere. 
Although the ripples found in the paper model of Ref. [5] are recovered, the overall structure 
with periodic boundary conditions is much flatter. 

These simple examples teach us that simulations of extensive grain boundary networks 
would be quite challenging. Finite size effects would make it difficult to discern the "real" 
effect, due to the breakdown of periodicity at a grain boundary vertex, from the spurious 
effect of freely terminating each of the boundaries at the edge of the sample. Note, however, 
that the leading contribution to the vertex stress tensor is the quadrupole moment. A bona 
fide localized defect possessing the same leading multipole stress field is a tightly bound 
dislocation pair. This is indeed amenable to computer simulation. 
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The buckling of disclinations and dislocations in a tethered membrane has been inves- 
tigated numerically by several authors [2,13]. Both types of defects induce large deflections 
in the membrane upon buckling. Figure 6 shows the buckling of a tightly bound dislocation 
pair. One can observe large out-of- plane displacements. 

An interesting question to address concerns the shape of a membrane in the presence of 
more than one grain boundary vertex, and the interaction which arises. An argument can 
be made that the buckling of two defects in opposite directions is favored by the bending 
energy (the stretching energy is less important in the buckled phase) [4]. Therefore, one 
may expect that point defects have an "antiferromagnetic" interaction. The simplest 
system where this prediction can be tested is a pair of impurities consisting of atoms of a 
bigger size, each of which produces a spherically symmetric compression of the neighboring 
lattice in the flat phase, with stresses that decay as 1/r^. For computational purposes, these 
defects were embedded near the center of an approximately circular membrane (radius 
R = 20 lattice spacings, and Kq/k 92, well above the buckling transition) and the 
energy was minimized by constraining the defects to buckle either on the same side or on 
opposite sides of the flat reference configuration. Bonds radiating from the impurity atoms 
had equilibrium length 1.5 times the bonds in the matrix. The energy difference between 
the two cases is taken as an estimate of the coupling (the exchange J of an Ising model). 
We found that J is positive for separations below three lattice constants, and negative 
above that (see Fig. 7). Thus, nearest neighbor and next-nearest neighbor impurities 
have a ferromagnetic coupling, which favors buckling in the same direction, while defects 
farther apart have an antiferromagnetic interaction, consistent with the argument above. 
Note the similarity with the frustrated interactions which arise in spin glasses [8]. In 
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real grain boundary networks, however, we have the additional possibility of a competing 
ferromagnetic interaction between nodes connected by a buckled grain boundary. 

One final feature may be important in the case of grain boundary vertices, namely, 
significant anisotropy of the exchange coupling. While detailed numerical study of the 
coupling in this more complicated case is beyond the scope of the present work, we present 
the result of a simulation of pairs of interacting defects, to show that the basic conclusion, 
i.e., an antiferromagnetic interaction between defects, is still valid. Figure 8 displays 
the topology of a membrane containing two "nodes," each modeled as a tightly bound 
dislocation pair. Each pair can be viewed as a collapsed (i.e., anisotropic) vacancy in 
a triangular lattice. The ratio of the elastic constants is Kob'^ / k pa 92, well above the 
buckling transition of an isolated pair, which we have found to occur for KqIP' / k 26. The 
relaxation was effected subject to the constraint that the 5-fold coordinated atoms marked 
by thicker circles in Fig. 8 buckle in the same (Fig. 9a) or in opposite (Fig. 9b) direction. 
The five-fold atoms nearest to these heavy circles buckle oppositely to their neighbors 
as in Fig. 6. The total energy is 9% smaller in Fig. 9b, indicating an antiferromagnetic 
interaction between these defects for this separation. 

V. DISCUSSION 

In conclusion, we hope to have clarified the way in which frustration arises in grain 
boundary networks. Particularly important are the nodes, which lead to long range strain 
fields and large out-of-plane displacements upon buckling. An infinite grain boundary 
by itself produces relatively small, exponentially screened defiections below the buckling 
transition. The stronger "antiferromagnetic" interaction between the buckled nodes across 
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a grain competes with a ferromagnetic nodal interaction which presumably arises when a 
buckled grain boundary connects two nearest neighbor nodes. 

In studies of spin glasses, one often considers simplified theoretical models, abstracted 
from the complexities of real quantum mechanical spins interacting via, say, a RKKY 
interaction [8] . Similar simplifications may be appropriate for grain boundary networks in 
membranes. A "direct" simulation of this problem might involve quenching a flat space 
liquid to produce grains, tethering neighboring particles to prevent further changes in 
topology, and then annealing and relaxing the resulting membrane in three dimensions. 
The free boundary conditions natural in this approach would lead to additional long range 
strain fields (and further warping of the surface) whenever a boundary terminates at the 
membrane edge. Sample sizes much larger than the grain would be required to accurately 
model the experiment of Ref. [3] , although even relatively small systems would be of some 
interest. One might also try to simulate grain boundary networks in membranes with a 
spherical topology, similar to the liposomes of Ref. [3] . 

In view of the complexity of direct simulations, it may be worth considering simplified 
models. One approach would be to deliberately insert vacancies and interstitials into an 
otherwise perfect membrane (to simulate nodes), and ignore the exponentially screened 
grain boundaries entirely. Such a model would be similar to the random impurity model 
discussed in Ref. [4]. For weak disorder, the equilibrium fiat phase [2] survives at finite 
temperatures (with microscopic buckling at small scales) , although there are indications of 
a spin glass phase when the disorder is strong [14]. Glassy behavior at low temperatures 
has been found in computer simulations [15]. One could abstract even further and simply 
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replace the buckled defects by Ising degrees of freedom with effective Hamiltonian 



where ai = ±1, the {fi} are random positions within the membrane, and J(r) has the form 
shown in Fig. 7. Simple power counting using Eq. (2) suggests that the antiferromagnetic 
tail of J(r) has the form J(r) ^ —n/r'^. Even the behavior of this simple two-dimensional 
Ising spin glass is, to the best of our knowledge, unknown. 
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Figure captions 

Fig. 1: The temperature behavior of the elastic constants of a membrane can cause finite 
energy defects, such as vacancies, interstitials, or grain boundaries, to buckle as the tem- 
perature is lowered. Buckling will occur at higher temperatures for defects characterized 
by a larger length scale, according to Eq. (1). Thus, low angle grain boundaries will buckle 
first upon cooling a polymerized grain boundary network 

Fig. 2: Closure failure for a small angle (9.4°) grain boundary indicating the presence of 
dislocations. 

Fig. 3: Symmetric high angle (21.8°) grain boundary in triangular lattice. The arrows 
display the points where the lattice is most strained. The pattern repeats itself indefinitely 
in the direction of the grain boundary. 

Fig. 4: Asymmetric high angle grain boundary in triangular lattice (see also Ref. [5]). For 
the grains to match, a uniform compressive load must be applied to the untilted boundary. 

Fig. 5: a: Buckling of a 21.8° symmetric grain boundary (the configuration in fiat space 
is that of Fig. 3a). Periodic boundary conditions have been applied in the direction of the 
grain boundary, b: Same as (a), except free boundary conditions have been applied in all 
directions, c: Buckling of asymmetric grain boundary (see Fig. 4 for the configuration in 
fiat space). Periodic boundary conditions have been applied in the direction of the grain 
boundary. 

Fig. 6: Buckled tightly bound pair of dislocations with opposite Burger's vectors (the 
configuration in fiat space is shown in inset). 
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Fig. 7: Energy splitting between "antiferromagnetic" and "ferromagnetic" configurations 
of two impurities in an otherwise ideal crystalline membrane. Each impurity is assigned a 
preferred bond length which is 50% larger than the lattice constant. An isolated impurity 
has a buckling transition for Kq/k, >8; here, we have taken Kq/k, = 92 (the lattice spacing 
is set to unity). When the impurities are farther apart than two lattice sites, buckling in 
opposite direction is energetically favored. The energy is given in units of the bending 
rigidity k. 

Fig. 8: Two tightly bound dislocation pairs in flat space. Explicitly marked are the 5-fold 
and 7-fold coordinated atoms. The unevenly spaced atoms reside in regions of high local 
strain. 

Fig. 9a,b: Relaxed configuration of a membrane comprising two tightly bound dislocation 
pairs (see Fig. 8 for the unrelaxed configuration). The relaxation was effected subject to 
the constraint that the 5-fold coordinated atoms (marked by a thicker circle in Fig. 8) 
buckle (a) in the same or (b) in opposite direction. The total energy is 9% smaller in the 
latter case, indicating an "antiferromagnetic" interaction between defects. 
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